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We give a Hamiltonian analysis of the asymptotically flat spherically symmetric 
system of gravity coupled to a scalar field. This 1+1 dimensional field theory may 
be viewed as the ’’standard model” for studying black hole physics. Our analysis 
is adapted to the flat slice Painleve-Gullstrand coordinates. We give a Hamiltonian 
action principle for this system, which yields an asymptotic mass formula. We then 
perform a time gauge fixing that gives a Hamiltonian as the integral of a local 
density. The Hamiltonian takes a relatively simple form compared to earlier work in 
Schwarzschild gauge, and therefore provides a setting amenable to full quantisation. 
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INTRODUCTION 


It may be argued that the central unsolved problem in black hole physics is a full quantiza¬ 
tion of the gravity-scalar held system in spherical symmetry. All the puzzles associated with 
black holes originated from studying this system classically and semi-classically m q 
The extra step to full quantization has never been attempted, mainly due to the intractabil¬ 
ity of the Hamiltonian system. If this could be accomplished, we would have a complete 
scenario for studying the formation and evaporation of Schwarzschild black holes in a fully 
quantum dynamical setting, and in the simplest ”no frills” context. 

Although there has been some progress in this area in string theory, it is restricted to 
highly special hlack holes in denved supecgnavhy models Q. In the loop quantum gtatnW 
program, the work on this problem has so far been restricted to ’’isolated” horizons p|. 
As the name suggests, these horizons are not appropriate for studying the full dynamics of 
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matter and gravity, where horizons can form and evolve in response to matter flows. 

The first attempt at starting a quantization of the gravity-scalar held model in four 
spacetime dimensions dates back to the seventies, when the Hamiltonian theory was worked 


0 , 


out in a parametrization adapted to Schwarzschild coordinates by Berger et. ah (BCNM) 
and later clarified by Unruh 2]- In this setting a reduced Hamiltonian for the scalar held was 
obtained in a particular gauge. This Hamiltonian was complicated enough that quantization 
was effectively untenable. (Unruh comments on this by saying ’’I present it here in the 
hope that someone else may be able to do something with it.” (2|). The prob 


since remained largely unaddressed, except for related work on shell collapse 




lid. 


em has 


12 1, and a 


geometrodynamical quantization of the (vacuum) Schwarzschild black hole 

In this paper we reanalyse the spherically symmetric gravity-scalar held system. The 
main new ingredient is the application of the Arnowitt-Deser-Misner (ADM) Hamiltonian 
ormulation to coordinates adapted to the hat slice Painleve-Gullstrand (PG) coordinates 


1 11 112L 11 Bl 114 nfl ]. with a time gauge fixing. The black hole metric in these coordinates is 
given by 


ds 2 = —dt 2 + \ dr + 


2 M 


dt + r^dQ- 


( 1 ) 


The spatial metric e a b given by the constant t slices is hat, and the extrinsic curvature of 
the slices is 

K n h \! ~ ( Gab ^ S a Sb ) , ( 2 ) 


where s a = x a /r for Cartesian coordinates x a 

The black hole mass information is contained only in the extrinsic curvature, which in 
the canonical ADM variables ( q a b , ir ab ) determines the momenta 7r ab conjugate to the spatial 
metric q a b. In this form, the mass formula is necessarily different from the ADM mass 
integral, since the spatial slices are hat. 

The PG coordinates motivate the following prescription for the falloff conditions for 
asymptotic flatness: 


q a b - e ab + +0(r 2 ) (3) 

7r ab = 9ab ^ + O (r -3//2-e ) , (4) 

where e > 0, f ab ,g ab ,h ab are symmetric tensors, TT ab = n ah /^/q, and q = det(g a fe). I 11 this 
definition it is manifest that the leading terms correspond to the black hole solution in PG 
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coordinates. The form of the next terms containing e is necessitated by the requirement of 
a well defined action principle. It is useful to compare these with the conditions motivated 
by the Schwarzschild coordinates, where the eading order terms in the metric and extrinsic 

nfl 

curvature are 1/r and 1/r 2 , respectively [17LIl8i]. 

These falloff conditions provide the starting point for our analysis of the Hamiltonian 
dynamics for spherically symmetric metrics of the form 


ds 2 = —f(r, t) 2 dt 2 + (dr + g(r, t)dt) 2 + r 2 dQ 2 , 


( 5 ) 


minimally coupled to a massless scalar field 0(r, t). 

In the next section, we give a parametrization for the ADM variables (q a b, n ab ) respecting 
these conditions to obtain a Hamiltonian theory. In Section III, we utilise a special time 
gauge fixing condition in which the reduced Hamiltonian takes a relatively simple form 
compared to the earlier works mentioned above. The setting therefore provides an arena in 
which a full quantization appears to be possible. The last section contains a brief comparison 
with the BCNM work, and a discussion of our approach to quantization. The present paper 
also provides the classical details which underlie the recent work on quantization by the 
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II. GRAVITY-SCALAR FIELD MODEL 


A well defined variational principle for Einstein’s equations coupled to matter, and sat¬ 
isfying a specified class of boundary conditions, may be obtained by starting with the bulk 
Einstein-Hilbert action, or its canonical bulk ADM form. The first requirement is that all 
the terms in the action are well-defined. The second requirement is that the variation of the 
action be of a form such that all variational derivatives with respect to the field variables 
are well defined. As noted by Regge and Teitelboim [l3, G|, this requires in general the 
addition of a surface term to the original action. 

The phase space of the model is defined by prescribing a form of the gravitational phase 
space variables q a b and ir ab , together with falloff conditions in r for these variables, and for 
the lapse and shift functions N and N a , such that the ADM 3+1 action for general relativity 
minimally coupled to a massless scalar field 

1 


S = 


167 tG 


d xdt 


n ab q a b + 


-NH - N a C a 


( 6 ) 
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is well defined. The constraints arising from varying the lapse and shift are 

n = 4= 

V 

-SvrG (— 

C a = D c n c a - 8nGP ri> d a (j) = 0, (8) 


- Vq r (q) 

+ + s/q^d^dtA = o (7) 


where n = n ab q a b■ The corresponding conditions for the matter fields (j) and P^ are deter¬ 
mined by the constraint equations. 

In this setting, the following parametrization for the 3-metric and conjugate momentum 
gives a reduction to spherical symmetry: 


Qab Sa,Sb H - 

■l ab 


R{r, tf 


(f+ft S a Sb) 


IT = 


Pa (jyt) 
2 A (r, t ) 


s a s b + 


r 2 P R (r,t) 
4 R(r, t ) 


(e ab - s a s b ) 


(9) 

( 10 ) 


Substituting these into the 3+1 ADM action shows that the pairs ( R , P R ) and (A, P\) 
are canonically conjugate variables. The reduced ADM 1+1 field theory action takes the 
form 


Sr 


1 

4 G 


dtdr (PrR + P\A + 


1 

4 G 


dtdr (NH + N r C r ) 


+ surface term, 


( 11 ) 


where we have performed the angular integral. The surface term is derived below. The 
Hamiltonian and diffeomorphism constraints are 


H 


C r 


1 

wx 

2 


I(f> A A) 2 -l(P A A)(f> R i?) 

[2RR"A - 2RR'A' - A 3 + A R' 2 ] 


n , r2 

2A R 2 2A 
PrR' ~ A P' K + 


U 2 


0. 


( 12 ) 

(13) 


These constraints are first class. The falloff conditions induced on the reduced variables by 
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(TO are 


R = r + 0(r^ 2 ~ e ), 

P R = Ar~ 1/2 /2 + 0(r- 1 ~ e ), 
A = l + C>(r- 3/2 - e ), 

P A = Ar 1 / 2 + 0(r~ e ) 

0 = Br~ 1/2 + 0(r~ 3/2 - e ), 

P,* = Cr 1 / 2 + 0(r- £ ), 


(14) 

(15) 

(16) 


where A, B,C are constants. This means that the asymptotic region is not dynamical (as 
it should be since it is flat). The constant A, which appears in the expressions for Pr and 
Pa above, will turn out to be captured by a surface integral (see below) and is proportional 
to the mass of the system. The above conditions together with the falloff conditions on the 
lapse and shift functions 

N r = Ar~ 1/2 + 0{r~ 1/2 ~ e ) 

N = 1 + 0(r ~ e ) (17) 


ensure that Sr is well-defined. More explicitly, they ensure that the symplectic structure is 
well-defined, which means that the integral of the terms PrR etc. converges. Furthermore, 
they guarantee that H vanishes to first order, and has a falloff beyond leading order, as 
required for the action to be well-defined. The same holds for the diffeomorphism constraint. 
Note that the factor 1/2 in the leading order term in Pr ensures that to this order both 
the diffeomorphism and Hamiltonian constraints vanish. This corresponds to the black hole 
solution in PG coordinates, which motivated our falloff conditions in the first place. Taken 
together these observations guarantee that the bulk action is well-defined. 

Consider now the variation of this action to see what surface terms need to be added to 


0 . 


make the variational principle well-defined. Following [17j, we compute the variation 5Sr, see 
what surface terms arise in it, and identify the terms that vanish due to the falloff conditions; 
the ones that do not must be compensated for by adding a surface term to the starting 
bulk action. In our case surface terms arise from those terms in the action that contain r 
derivatives. These are the Ricci scalar, matter density, and the radial diffeomorphism terms. 
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The variation is 


4 G 5Sr = I dtdr (terms giving eqns. of motion) 

NR 2 (j) r 


dt 

dt 

dt 


N r Pa, + 


A 


8(f> 


ANRR' 

N r A 5P A + -—- 8A 


NR 


A J A 
- I dt [N r P R } 5R, 


A 2 


AN R 1 4NRN 

+ 


A 2 


8R 


(18) 


where the last four terms are differences of surface integrals evaluated at r — 0 and r — oo. 
The variational principle is well-defined if each of these surface terms vanishes. 

At r = oo most of these vanish by virtue of the falloff conditions. This leaves two terms 
of order one. The one proportional to <50 can be eliminated by requiring this variation to 
vanish at infinity. The other, proportional to 6P\, is dealt with by adding a surface term at 
infinity to the original bulk action whose variation cancels the offending surface term from 
the variation. It is this term that captures the conserved asymptotic mass. 

At r = 0 we do not impose falloff conditions on the phase space variables because there 
is no physical guidance for this. We thus require the usual prescription that the variations 
of the configuration variables vanish there: 


8(f)\r=0 0, 

dR\r=0 = 0, 

5A| r=0 = 0. (19) 


This leaves only the term proportional to 5P\ at r = 0. To deal with this we require either 
the addition of a surface term with the opposite sign to the one at infinity, or the condition 
N r (r = 0) = 0. The former would subtract from the asymptotic mass. Therefore for the 
vacuum solutions (0 = 0 and = 0) in the gauge R = r A = 1, the sum of the surface 
terms at r = 0 and r — oo would cancel. This suggests that we make the latter choice. 

Based on these observations, functional differentiability of the action is guaranteed if we 
add the term 

- J dt(N r AP A ) 


r=oo 


( 20 ) 
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to the original bulk action. We can then derive the evolution equations: 


R 

Pr 


A 

Pa 


P+ 


-N ^ + N r R' 
4 R 


N 


V Pl A 


PrPa + Pi 


R 


J 2 


4P 3 4P 2 AP 3 A 

/4PAWV 4PW' /4P7VY 


A 


V A y 


V A 2 

+(WP R )' 

(ftA 2 - A rp r ) + (A rry 

N 1 Ja 


Pi 4 RR" „ 2(P') 2 


8R 2 
fARR’N 

~ v A 2 
+N r P' k 


+ 


A 2 


+ 


+ 2 + 


8PPW 


A 2 


A 3 


iV 


2A 2 R 2 


(P| + P 4 (0O 2 ) 


AT 


AP 2 


+ ivy 


AT 


yp 2 

T“ 


+ (iV r P 0 )'. 


( 21 ) 


( 22 ) 

(23) 


(24) 

(25) 

(26) 


Note that the surface term d25D is the mass formula for the flat slice parametrization. 
Substituting the asymptotic forms of the variables m and m gives 


JVAPa = A 2 + O(r-), 


(27) 


which shows that the parameter A contains conserved mass information. The relation be¬ 
tween A and the conventional mass parameter M in (HD is obtained by comparing the 7r ab 
obtained from @ with our parametrization m ■ This gives A = 4v/2 M. We note also 
that our falloff conditions are preserved under this evolution, which ensures their consis¬ 
tency. This is easily seen by computing the left and right hand sides of the above evolution 
equations in the asymptotic regime. 


III. TIME GAUGE FIXING 

We now gauge fix the theory defined above with the condition A = 1. It is second class 
with the Hamiltonian constraint. Demanding that it be preserved in time implies from (ESI) 
the relation 


N(P a - RP r ) = -AR 2 {N r )' 


( 28 ) 
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between the lapse N and the shift N r . As a result of the time gauge fixing the Hamiltonian 
constraint m must be imposed strongly. This gives a quadratic equation for Pa in terms of 
the remaining variables. A comparison with the vacuum solution in fully gauge-fixed form 
(ie. with the coordinate fixing condition R = r), uniqely selects the positive root. This gives 


P A = P R P + V(/V?) 2 -A', 


(29) 


where 

X = 16P 2 (2PP" - 1 + R' 2 ) + 16 R 2 H (t> (30) 


and 

h * = w + t ^ (31) 

The positivity of the argument of the square root follows from the dominant energy condition 
for the massless scalar field. 

The solution for the lapse function now reads 


AR 2 (N r )' 

v/(P s P) 2 -.Y' 


(32) 


The reduced Hamiltonian equations for the remaining canonical variables (P, Pr) and 
( 0 , Pfj,) are obtained by substituting the gauge condition A = 1, the corresponding solution 
© of the Hamiltonian constraint, and the lapse equation (13211 into the unfixed evolution 
equations (EH)-©, and into the radial cliffeomorphism constraint. The gauge fixed equa¬ 
tions are 


R = ~ (PrP+VCPrA) 2 - v) + N T R' 
Pr = (PrR 2 + PrPV(PrP ) 2 - V - A') 

+N (A _ Rr:p \ + 47J'AT' _ (4P.Y)' 
+(N r p R y 
i = N A + JVV 
p» = (iv^'p 2 )' + (N r p t y, 


(33) 


(34) 

(35) 

(36) 


where it is understood that N is given by m- The remaining radial diffeomorphism con¬ 


straint is 


P\ + PrR' + Pfytf — 0 , 


( 37 ) 









9 


where P\ is given by (BUD- 

All these equations can be obtained from the gauge fixed reduced action 


cG _ 
D R — 


dtdr 


+ p r r - iv r (p; + p r r! + p$') 


dtp A (N r y. 


(38) 


The surface term may be combined with the bulk term to write the action in a form from 
which one can read off the gauge fixed Hamiltonian 


Hr = 


[(N r )'P A + N r (P R R' + Pt/tcf)')} dr 
(N r y {RP r + y/{P R R)2-X) dr 

/»oo 

1 N r (P R R' + P^') dr, 


(39) 


IV. DISCUSSION 


The Hamiltonian © is a simpler expression than that obtained from the full time and 
coordinate gauge fixing given in |2j, where Schwarzschild gauge is used. That Hamiltonian 
is 

H s = I dr + r 2 (j) ,2 \ exp ( j S (j> (r')dr 1 ') (40) 

where 

P? r<ifl 

(41) 


„,, n , ^ 

S *(r) = -^ + — 


It is apparent that the square root in m is easier to handle than the non-locality manifest 
in this formula. 

The non-local term can be traced back to the fact that (TO is arrived at after both a time 
and radial coordinate gauge fixing. This may be seen in our formulation as well: Consider 
the radial gauge fixing R(r, t ) = r. This leads to (i) fixing of the shift function N r by the 
condition 

rPrt 


r = 0 = (N r )'r 1 + 


+ N r 


(42) 


y(r P R f-X 

and (ii) fixing of Pr by strong imposition of the remaining diffeomorphism constraint (HUD- 
The solution of (H2I) for N r is 

y/(rP R )*-X 


N r = exp 


fV(rP R ) 2 -X + f*P R 


dr 


(43) 
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which leads directly to a non-local term when substituted back into (HP- Thus we learn 
that the complicated form of the fully gauge fixed Hamiltonian is traceable solely to the 


le time gauge fixing. A partial gauge fixing in the diagonal 
would also lead to a simpler, local form of the reduced 


coordinate gauge fixing, but nott 
parametrization of the metric 
Hamiltonian. 

For this reason, we propose that a quantization of this system be carried out in a par¬ 
tially gauge fixed setting, with only the time gauge fixed. Retaining a first class radial 
diffeomorphism constraint presents issues that are relatively easier to deal with; we have in 


fact already proposed Q] a kinematical quantization of this system in which the Hilbert 
space carries a representation of finite diffeomorphisms. 

In summary, we have presented here in detail the canonical formalism for spherically 
symmetric gravity coupled to a scalar held adapted to the hat slice foliation. We have 
shown that all consistency conditions are satished. These include functional differentiability, 
surface terms, and the preservation under evolution of the falloff conditions at infinity. 

Having obtained the reduced Hamiltonian m, the next step in our program is to con¬ 


struct t 
appear 


re corresponding operator using the techniques presented 


in jwj. 


This work is to 


|2l|. 
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